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Abstract 

Background: Expansion of aquaculture requires alternative feeds and breeding strategies to reduce dependency 
on fish oil (FO) and better utilization of dietary vegetable oil (VO). Despite the central role of intestine in 
maintaining body homeostasis and health, its molecular response to replacement of dietary FO by VO has been 
little investigated. This study employed transcriptomic and proteomic analyses to study effects of dietary VO in two 
family groups of Atlantic salmon selected for flesh lipid content, 'Lean' or 'Fat'. 

Results: Metabolism, particularly of lipid and energy, was the functional category most affected by diet. Important 
effects were also measured in ribosomal proteins and signalling. The long-chain polyunsaturated fatty acid 
(LC-PUFA) biosynthesis pathway, assessed by fatty acid composition and gene expression, was influenced by 
genotype. Intestinal tissue contents of docosahexaenoic acid were equivalent in Lean salmon fed either a FO or VO 
diet and expression of LC-PUFA biosynthesis genes was up-regulated in VO-fed fish in Fat salmon. Dietary VO 
increased lipogenesis in Lean fish, assessed by expression of FAS, while no effect was observed on (3-oxidation 
although transcripts of the mitochondrial respiratory chain were down-regulated, suggesting less active energetic 
metabolism in fish fed VO. In contrast, dietary VO up-regulated genes and proteins involved in detoxification, 
antioxidant defence and apoptosis, which could be associated with higher levels of polycyclic aromatic 
hydrocarbons in this diet. Regarding genotype, the following pathways were identified as being differentially 
affected: proteasomal proteolysis, response to oxidative and cellular stress (xenobiotic and oxidant metabolism and 
heat shock proteins), apoptosis and structural proteins particularly associated with tissue contractile properties. 
Genotype effects were accentuated by dietary VO. 

Conclusions: Intestinal metabolism was affected by diet and genotype. Lean fish may have higher responsiveness 
to low dietary n-3 LC-PUFA, up-regulating the biosynthetic pathway when fed dietary VO. As global aquaculture 
searches for alternative oils for feeds, this study alerts to the potential of VO introducing contaminants and 
demonstrates the detoxifying role of intestine. Finally, data indicate genotype-specific responses in the intestinal 
transcriptome and proteome to dietary VO, including possibly structural properties of the intestinal layer and 
defence against cellular stress, with Lean fish being more susceptible to diet-induced oxidative stress. 
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Background 

Fish are important components of the human diet, being 
highly nutritious and valued as the main source of n-3 
long-chain polyunsaturated fatty acids (LC-PUFA). 
These essential fatty acids, mainly eicosapentaenoic acid 
(EPA) and docosahexaenoic acid (DHA), have well- 
known health-promoting properties, including protec- 
tion against a range of cardiovascular and inflammatory 
diseases, and neurological disorders [1]. With population 
growth and increasing awareness of the importance of 
fish consumption as part of a healthy diet, worldwide de- 
mand for seafood continues to grow. However, as trad- 
itional fisheries are largely in decline, aquaculture must 
meet this demand. Aquaculture is the fastest-growing 
food production sector with an average annual growth 
rate of 6.6%, accounting for 46% of total fish supply [2]. 
In the European and American continents, aquaculture 
production is largely dominated by salmonid species, 
mainly Atlantic salmon, and feeds for such carnivorous 
species have traditionally relied on fishmeal (FM) and 
fish oil (FO) from wild stocks. Recent estimates indi- 
cated that 88.5% of global production of FO was used by 
the aquaculture sector, with salmonid culture taking the 
largest share (56%) [3]. With ever increasing demands 
for aquafeeds and reduction in fisheries landings, the 
availability of FO and FM seriously limits the growth of 
aquaculture production and there is an urgent need to 
find more sustainable alternatives. 

Vegetable oils (VO) can replace FO in salmon feeds 
without compromising fish growth or condition al- 
though, at high levels of replacement, tissue levels of n-3 
LC-PUFA are significantly reduced [4,5]. The effects of 
FO replacement by VO are becoming well characterized 
in the hepatic transcriptome of salmonids [6-10], and 
other species [11]. However, studies on intestinal tran- 
scriptome are few and restricted to effects of replace- 
ment of FM by plant proteins, particularly soybean meal 
[12,13], given its potential to cause enteritis. It is now 
clear that intestine in salmonids is not simply a site for 
reacylation and packaging of dietary lipids but it also has 
important roles in fatty acid metabolism, including LC- 
PUFA biosynthesis [14,15]. Furthermore, dietary VO can 
induce major histological changes in fish enter ocytes, 
originating mostly from supranuclear lipid droplet for- 
mation, possibly due to altered reacylation mechanisms 
and decreased phospholipid synthesis [16-18]. In some 
cases, these accumulations were large enough to be 
deemed pathological [19,20]. A recent study investigat- 
ing effects of dietary FO replacement by VO on intes- 
tinal transcriptome in Atlantic cod indicated potential 
effects on lipid absorption and transport and suggested 
morphological and structural changes to the intestinal 
muscle layer [21]. Furthermore, both this and a previous 
study on Atlantic salmon showed significant effects on 



expression of genes involved in cell proliferation and 
apoptosis [22]. Therefore, there is indication that intes- 
tine may be affected by changes in lipid components of 
feed formulations. Given its crucial roles in nutrient ab- 
sorption, protection against the entry of pathogens, and 
immune function [23], further attention is warranted 
and impacts of FO replacement require investigation in 
intestine, particularly in salmon where important 
changes in diet formulation are already being applied. 
This study is a large-scale analysis of the effects of re- 
placement of dietary FO by VO on the transcriptome 
and proteome of Atlantic salmon intestine. Furthermore, 
given recent interest in evaluating genetic selection as a 
feasible strategy, in conjunction with changes in com- 
mercial feed formulation, to meet worldwide demand for 
farmed fish without compromising animal welfare or nu- 
tritional value [9,10,24], two groups of Atlantic salmon 
families, Lean and Fat, were studied to examine the po- 
tential effects of genetic background. This experiment 
was performed in parallel with another microarray study 
looking at effects in the hepatic transcriptome, analysing 
samples from the same individuals [9,10], enabling a glo- 
bal and comprehensive assessment of the physiological 
and molecular effects of FO replacement by VO in At- 
lantic salmon, including potential interactions with 
genotype. 

Results 

Microarray analysis 

Two way-ANOVA of the cDNA array dataset returned 
1409, 1626 and 862 significant genes for the factors diet, 
genotype and diet x genotype interaction, respectively. 
Detailed analysis was restricted to the top 100 most sig- 
nificant features, which were categorised according to 
biological function, based on mammalian homolog 
genes. Metabolism, particularly of lipid and energy, was 
the functional category most affected by diet accounting 
for 39-41% of the top 100 annotated genes (Table 1), 
and showing highest diet x genotype interaction (Add- 
itional file 1). Diet also impacted translation (17%; 
mostly ribosomal proteins) and signalling (17%). In con- 
trast, genotype affected less markedly metabolism (21%), 
whereas structural proteins (25%; mostly different types 
of collagen) and proteins involved in the regulation of 
transcription (21%) predominated (Table 2). 

Gene Ontology (GO) enrichment analysis was per- 
formed on the complete significant lists, enabling identi- 
fication of GO terms significantly enriched in the input 
entity list, in comparison to the whole array, providing 
clues as to which biological processes might be particu- 
larly altered in the experimental conditions being com- 
pared. It revealed no significant enrichment of GO 
terms in the genotype list, while 20 and 7 GO terms 
were significantly enriched in the diet and interaction 
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Table 1 Intestine transcripts corresponding to the top 100 most significant features exhibiting differential expression 
between diets 



Accession no 



Gene 



VO/FO 
Lean Fat 



p- value 



Metabolism (39%) 

Lipid metabolism (9%) 
DW590668 



CK884623 
CK885725 
AF478472 



Fatty aldehyde dehydrogenase 
Ethanolamine kinase 1 

Epidermis-type arachidonate lipoxygenase 3 
Protein containing a beta-ketoacyl synthase domain 
Delta-5 fatty acyl desaturase 



Energy metobolism/generotion of precursor metabolites (11%) 

CK885194 Cytochrome c oxidase subunit 5A 

bra_snb_03E12 NADH dehydrogenase subunit 1 

EG649067 Ubiquinol-cytochrome c reductase core protein 1 

CK884638 Uncoupling protein 

CK890974 Mitochondrial calcium-dependent solute carrier SCaMC-2 

EG648806 NADH dehydrogenase (ubiquinone) 1 beta subcomplex 9 

Protein and amino acid metabolism (9%) 

CK877521 Kunitz-type protease inhibitor 2 precursor 

DW590246 Betaine aldehyde dehydrogenase 

CK885180 Fumarylacetoacetate hydrolase domain-containing protein 2B 

CK897653 Ubiquitin-conjugating enzyme E2 A 

CK885145 Lysosomal protective protein precursor (Cathepsin A) 

Carbohydrate metabolism (5%) 

EG649263 Endosulfine alpha 

CK885545 Glucose transporter type 8 

Xenobiotic and oxidant metabolism (5%) 

CK885392 Catalase 

C0472279 Cytochrome P450 1 a 

CK898307 Selenoprotein H-like protein 

Transport/ intracellular trafficking (5%) 

EG647578 Calcium-activated potassium channel subunit alpha 1 

EG648841 Sodium-coupled neutral amino acid transporter 2 

EG648303 Mitochondrial import receptor subunit TOM20 homolog 

Regulation of transcription (13%) 
CK898304 
CK890154 
EG647831 
AM041851 
CK885179 
CK884728 
CK879187 
Translation (17%) 
CK878380 
DW591873 



C-Maf 

Butyrate response factor 2 
Proliferation-associated 2G4 b 
Enhancer trap locus homolog 1 
Steroid receptor RNA activator 1 
Butyrate response factor 1 
Enhancer of polycomb homolog 1 

40S ribosomal protein S18 
60S acidic ribosomal protein 



1.1 
1.2 
1.4 
1.4 
1.2 

1.2 
1.3 
1.2 
1.3 
1.5 
1.1 

1.2 
1.4 
1.4 
1.2 
1.1 

1.3 
1.0 



1.2 
1.3 

1.3 
1.2 
1.1 

1.6 
1.6 
1.3 
1.1 
1.3 
1.1 
1.4 

1.1 
1.3 



■ 1.2 
1.2 
1.1 
1.6 
2.5 

• 1.1 

■ 1.2 

■ 1.1 
-1.3 

1.0 

■ 1.1 

1.1 

■ 1.0 

■ 1.3 
1.2 
1.1 

■ 1.0 
1.5 

1.1 
1.4 
1.1 

1.2 

■ 1.2 

■ 1.1 

1.6 
1.2 

■ 1.2 

■ 1.2 

■ 1.3 
1.4 
1.1 

■ 1.2 

■ 1.2 



0.0003 
0.0008 
0.0010 
0.0017 
0.0019 

0.0000 
0.0000 
0.0004 
0.0019 
0.0021 
0.0022 

0.0002 
0.0005 
0.0007 
0.0012 
0.0014 

0.0005 
0.0013 

0.0000 
0.0004 
0.0013 

0.0005 
0.0007 
0.001 1 

0.0000 
0.0001 
0.0003 
0.0004 
0.0007 
0.001 1 
0.0020 

0.0003 
0.0006 
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Table 1 Intestine transcripts corresponding to the top 100 most significant features exhibiting differential expression 
between diets (Continued) 



EG648644 


60S ribosomal protein L29 


- 1.1 


- 1.2 


0.001 1 


DW590580 


40S ribosomal protein S5 


- 1.2 


- 1.2 


0.001 1 


CK880046 


40S ribosomal protein S28 


- 1.1 


- 1.2 


0.0017 


EG649167 


60S ribosomal protein L37a 


- 1.1 


- 1.2 


0.0017 


EG648512 


60S ribosomal protein L41 


- 1.1 


- 1.2 


0.0017 


CK886042 


60S ribosomal protein L31 


- 1.1 


- 1.3 


0.0021 


CK895350 


60S ribosomal protein L36 


- 1.1 


- 1.1 


0.0021 


Signalling/Signal transduction (17%) 








CK874863 


Proto-oncogene serine/threonine-protein kinase 


1.1 


1.1 


0.0002 


CK881266 


Mitogen-activated protein kinase kinase kinase kinase 1 


1.2 


1.2 


0.0002 


CK878590 


Serine/threonine protein kinase (Tribbles homolog 2) 


- 1.1 


- 1.3 


0.0004 


AJ425351 


Tensin-like CI domain-containing 


- 2.0 


- 1.2 


0.0006 


CK885079 


Sonic hedgehog-like protein 


1.3 


1.0 


0.0008 


CN181161 


TGF-beta type II receptor 


1.1 


1.1 


0.0009 


EG649257 


Serine/threonine-protein kinase VRK3 


- 1.1 


- 1.2 


0.0009 


CK886267 


Guanylate cyclase 2C 


1.1 


1.1 


0.001 1 


CK884467 


Serine/threonine-protein phosphatase 2A regulatory subunit B 


- 1.4 


- 1.3 


0.0020 


Structural proteins (7%) 










mus_snm_13F09 


Myosin heavy chain 


- 1.2 


- 1.1 


0.0002 


CK8851 18 


Cardiac tropomyosin 


1.0 


1.1 


0.0013 


CK883354 


Troponin C 


1.2 


1.2 


0.0017 


EG647866 


Tubulin alpha chain 


- 1.2 


- 1.2 


0.0020 


Immune response (2%) 










BM414177 


Liver-expressed antimicrobial peptide 2 


1.5 


1.4 


0.0022 


Miscellaneous 










CK884634 


Ezrin-radixin-moesin binding phosphoprotein 50 


1.1 


1.1 


0.0003 


EG648821 


Transposase 


- 1.1 


- 1.2 


0.0004 


CK883634 


Poly [ADP-ribose] polymerase 12 


1.3 


1.1 


0.0006 


CK885090 


Claudin 15 


1.2 


1.1 


0.0008 


DW589291 


Transposase 


- 1.3 


- 1.3 


0.0009 


CK878888 


Splicing factor, arginine/serine-rich 3 


- 1.3 


- 1.2 


0.0013 


CK885098 


RNA-binding protein with multiple splicing 2 


- 1.4 


- 1.2 


0.0013 


CK891747 


Occludin 


- 1.2 


- 1.2 


0.0014 


CK888908 


All-trans-retinol 13,14-reductase precursor 


1.1 


1.2 


0.0017 


CK880622 


Cyclin G1 


- 1.1 


- 1.2 


0.0019 


CK885604 


Methylenetetrahydrofolate dehydrogenase (NADP+ dependent) 1 -like 


- 1.6 


- 1.3 


0.0020 


CK891699 


Barrier-to-autointeg ration factor 


- 1.1 


- 1.2 


0.0020 


CK879123 


Hemoglobin subunit beta 


1.5 


1.0 


0.0022 



Annotated features (67% of all clones) are arranged by categories of biological function (with percentages of distribution of genes within each category being 
indicated) and, within these, by decreasing significance (assessed by two-way ANOVA). Also indicated are the GenBank accession numbers for each clone (or, 
when not available, the probe number) and the expression ratio between fish fed VO and those fed FO, for each genotype (Lean and Fat). 



lists, respectively. GO terms enriched in the diet list expression, macromolecule and biopolymer biosynthetic 

included structural constituents of ribosome, structural process and other related terms. This was explained by 

molecule activity, cytosolic ribosome, cytosol, ribosomal the large number of ribosomal proteins, components of 

subunit, translation, cellular biosynthetic process, gene both the 40S and 60S subunits, which were down- 
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Table 2 Intestine transcripts corresponding to the top 100 most significant features exhibiting differential expression 
between family groups 

Gene Lean/Fat p- value 



Accession no 



Lean/Fat 
FO VO 



Metabolism (21%) 

Lipid metabolism (5%) 

CK874747 Acylglycerol kinase 

BI468033 ATP-binding cassette sub-family A member 1 

Energy metobolism/generotion of precursor metabolites (2%) 

CK890974 Mitochondrial calcium-dependent solute carrier SCaMC-2 

Protein and amino acid metabolism (12%) 

EG647485 Proteasome subunit beta type-8 precursor 

DW591 059 Heat shock protein Hsp-1 6 

AJ425553 26S proteasome non-ATPase regulatory subunit 9 

DW592187 26 proteasome complex subunit DSS1 

CK8781 09 Heat shock protein 60 kDa 

Xenobiotic and oxidant metabolism (2%) 

C0472279 Cytochrome P450 1A 

Transport/ intracellular trafficking (5%) 

CK886307 Peptide transporter PEPT2 

DW588251 Post-GPI attachment to proteins factor 2 

Regulation of transcription (21%) 

CK877881 Forkhead box Q1 

CK884548 Proline-rich nuclear receptor coactivator 2 

CK885013 Forkhead box protein F1 

CK891975 Y-box binding protein-2 

DW589341 Signal transducer/activator of transcription Statl 

CN 181 282 Bromodomain containing protein 3 

CK884407 C-Myc-binding protein 

CK895950 Transcription factor CP2-like 

EG648473 BTB/POZ domain-containing protein 

Translation (7%) 

CN 181 322 40S ribosomal protein S4 

CK885979 Ribonuclease UK1 14 

EG648042 Eukaryotic translation initiation factor 3, subunit J 



Signalling/Signal transduction (16%) 



C0471793 

CK885560 

CK885079 

EG648038 

DW589610 

CK897590 

EG647545 

Structural proteins (25%) 

CK882401 
CK889952 



Rho GDP-dissociation inhibitor 1 

Protein tyrosine phosphatase non-receptor type 9 

Sonic hedgehog-like protein 

Kalirin, RhoGEF kinase 

Chimaerin 1 

Growth factor receptor-bound protein 2 
Protein tyrosine phosphatase receptor type D 

Troponin C 
Transgelin 2 



- 1.0 

- 1.1 

1.1 

- 1.3 
1.2 

- 1.2 
1.5 
1.2 

1.2 

- 1.4 
1.1 

- 1.4 

- 1.1 
1.1 

- 1.1 

- 1.2 

- 1.1 

- 1.3 

- 1.5 

- 1.1 

1.1 

- 1.1 
1.1 

1.3 
1.3 

- 1.3 
1.2 
1.2 
1.1 
1.0 

- 1.1 
1.5 



- 1.1 

- 1.6 

1.6 

- 1.8 
1.4 

- 1.1 
1.3 
1.4 

1.5 

- 1.3 
1.3 

- 1.7 

- 1.2 
1.5 

- 1.1 

- 1.2 

- 1.1 

- 1.1 

- 1.2 

- 1.4 

1.1 

- 1.3 
1.1 

1.0 
1.5 

- 1.0 
1.2 
1.2 
1.1 
1.2 

- 1.2 
1.7 



0.0003 
0.0004 

0.0000 

0.0000 
0.0001 
0.0002 
0.0008 
0.0009 

0.0003 

0.0006 
0.0007 

0.0000 
0.0001 
0.0003 
0.0003 
0.0004 
0.0007 
0.0007 
0.0007 
0.0009 

0.0002 
0.0005 
0.0009 

0.0000 
0.0004 
0.0004 
0.0005 
0.0005 
0.0009 
0.0009 

0.0000 
0.0000 
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Table 2 Intestine transcripts corresponding to the top 100 most significant features exhibiting differential expression 
between family groups (Continued) 



EG648135 


Osteonectin 


1.3 


1.5 


0.0000 


CK899058 


Collagen alpha 3 type VI 


1.5 


1.5 


0.0002 


CK879405 


Collagen alpha 2 type VI 


1.3 


1.5 


0.0002 


EG649361 


Collagen alpha 2 type I 


1.5 


1.7 


0.0005 


EG649013 


Collagen alpha 3 type I 


1.4 


1.5 


0.0005 


mus_mfo_15A10 


Troponin-I isoform 3 


1.4 


1.6 


0.0007 


CK876833 


Keratocan 


- 1.0 


1.6 


0.0007 


DW590534 


Collagen alpha 2 type I 


1.5 


1.6 


0.0010 


CK892271 


Collagen alpha 2 type V 


1.2 


1.4 


0.00010 


CK873441 


Collagen alphal type VI 


1.2 


1.4 


0.0010 


Immune response (5%) 










AJ425527 


T-cell immunoglobulin and mucin domain-containing protein 4 


1.1 


1.1 


0.0001 


EG649194 


Mannose-binding protein C 


1.4 


1.4 


0.0003 


Miscellaneous 










AM041770 


Preimplantation protein 3 


- 1.3 


- 1.2 


0.0000 


EG649106 


Ornithine decarboxylase antizyme 2 


1.4 


1.1 


0.0002 


CK884355 


Guanylin precursor 


1.2 


1.4 


0.0002 


CK879185 


Aminolevulinate, delta-, synthetase 1 


-2.0 


- 1.1 


0.0002 


CK885814 


Similar to p53-associated parkin-like cytoplasmic protein 


1.2 


1.1 


0.0002 


CK886214 


Septin-7 


1.1 


1.1 


0.0003 


BM414079 


Tumor protein D52-like 2 


- 1.1 


- 1.3 


0.0003 


BM413818 


Splicing factor, proline- and glutamine-rich 


1.2 


1.2 


0.0003 


CK874550 


Follistatin-related protein 1 


1.1 


1.3 


0.0004 


DW589104 


Nucleoside diphosphate kinase 6 


1.1 


1.1 


0.0005 


DW592207 


Deoxyhypusine synthase 


1.5 


1.6 


0.0007 


CK892359 


HCLS1 -associated protein X-1 


- 1.2 


- 1.1 


0.0007 


BI468158 


Nuclear protein 1 


- 1.0 


- 1.6 


0.0009 


AJ425003 


Lamina-associated polypeptide 2 


1.2 


1.1 


0.0010 


CK899283 


Small EDRK-rich factor 2 


1.3 


1.5 


0.0010 


CK885090 


Claudin 15 


- 1.2 


- 1.1 


0.001 1 



Annotated features (60% of all clones) are arranged by categories of biological function (with percentages of distribution of genes within each category being 
indicated) and, within these, by decreasing significance (assessed by two-way ANOVA). Also indicated are the GenBank accession numbers for each clone (or, 
when not available, the probe number) and the expression ratio between Lean and Fat fish fed either FO or VO. 



regulated by dietary VO (Table 1). In contrast, several 
A6-desaturase (A6fad) clones showing a diet x genotype 
interaction (see Additional file 1) caused a significant en- 
richment of the GO terms oxidoreductase activity, 
stearoyl-CoA 9-desaturase activity, unsaturated fatty acid 
biosynthetic activity/metabolic processes and very long 
chain fatty acid (VLCFA) biosynthetic activity/metabolic 
processes. 

RT-qPCR analysis of gene expression 

The expression of several genes significantly affected or 
related to processes affected by the two factors in the 
microarray analysis was determined by RT-qPCR 
(Table 3). For diet, a reasonably good match was found 



for A5 fatty acyl desaturase (A5fad), NADH dehydrogen- 
ase subunit 1 (NADH1), proliferation-associated 2G4b 
(PA2G4), 60S acidic ribosomal protein (RP60S), prolifer- 
ating cell nuclear antigen (PCNA; also known as cyclin) 
and cytochrome P450 1A (CYP1A), particularly in the 
Fat group where fold-changes were generally more pro- 
nounced and significant. No change in expression of un- 
coupling protein 2 (UCP2) with diet was measured 
while, for myosin heavy chain (MYO) and methylenete- 
trahydrofolate dehydrogenase (NADP+ dependent) 1- 
like (MTHFD1), RT-qPCR indicated a change opposite 
to that suggested by microarray. Regarding genotype, a 
good match was obtained for CYP1A, proteasome sub- 
unit beta type-8 precursor (PSMB8) and alpha 2 type I 
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Table 3 Relative expression of genes assayed by RT-qPCR 



Genes VO/FO Lean/Fat 





Lean 






Fat 




FO 


VO 






Ratio 


p value 


Ratio 


p value 


Ratio 


p value 


Ratio 


p valu 


Microarray validation 


















A5fad (D) 


1.72 


ns 


4.34 


0.005 


2.49 


ns 


-1.01 


ns 


NADH1 (D) 


-1.03 


ns 


-1.28 


0.05 


1.00 


ns 


1.24 


ns 


UCP2 (D) 


-1.10 


ns 


1.08 


ns 


-1.10 


ns 


-1.28 


0.003 


MTHFD1 (D) 


1.20 


ns 


-1.09 


ns 


-1.08 


ns 


1.22 


ns 


PA2G4 (D) 


-1.27 


ns 


-1.23 


0.032 


1.11 


ns 


1.1 


ns 


RP60S (D) 


-1.02 


ns 


-1.69 


0.023 


1.13 


ns 


1.9 


0.006 


MYO (D) 


1.77 


0.034 


1.22 


ns 


-1.35 


ns 


1.08 


ns 


PCNA (D) 


-1.14 


ns 


-1.27 


ns 


1.08 


ns 


1.20 


ns 


CYP1A (D+G) 


1.71 


0.004 


1.81 


0.023 


-2.22 


0.003 


-2.38 


0.001 


ABCA1 (G) 


1.24 


ns 


1.17 


ns 


1.21 


ns 


1.28 


0.015 


PSMB8 (G) 


-100.00 


ns 


1.87 


0.033 


-16.67 


0.026 


< -1000 


0.000 


TAGLN (G) 


1.23 


ns 


-1.02 


ns 


-1.03 


ns 


1.22 


ns 


C0L1A2 (G) 


1.31 


ns 


-1.06 


ns 


1.47 


ns 


2.06 


0.001 


LC-PUFA biosynthesis and regulation 














A6fad_a 


1.11 


ns 


2.96 


0.008 


1.76 


ns 


-1.52 


ns 


elovl5a 


1.08 


ns 


1.2 


ns 


1.03 


ns 


-1.08 


ns 


elovl5b 


-1.19 


ns 


1.92 


0.031 


1.54 


ns 


-1.47 


ns 


elovl2 


1.2 


ns 


14.74 


0.003 


4.89 


ns 


-2.50 


ns 


Fatty acid synthesis 


















FAS 


1.49 


0.010 


1.20 


ns 


1.05 


ns 


1.31 


ns 


Fatty acid oxidation 


















ACO 


1.18 


ns 


1.28 


ns 


1.19 


ns 


1.09 


ns 


CPT1 


1.23 


ns 


-1.09 


ns 


-1.16 


ns 


1.15 


ns 


Regulation of lipid metabolism 
















PPARa 


1.14 


ns 


1.38 


ns 


1.26 


ns 


1.04 


ns 


PPARP 


1.03 


ns 


1.03 


ns 


-1.02 


ns 


-1.02 


ns 


PPARy 


1.52 


ns 


1.60 


0.036 


1.21 


ns 


1.15 


ns 


Xenobiotic and oxidative stress 
















GST 


-1.30 


ns 


1.07 


ns 


1.21 


ns 


-1.14 


ns 


MT-A 


-1.82 


o.05 r 


-1.05 


ns 


-1.25 


ns 


-2.13 


0.007 


CAT 


1.57 


0.036 


1.26 


ns 


-1.52 


ns 


-1.22 


ns 


SOD 


-1.20 


ns 


-1.04 


ns 


1.21 


ns 


1.04 


ns 


Apoptosis 


















CASP3B 


1.22 


0.045 


1.22 


ns 


-1.04 


ns 


-1.04 


ns 


CASP6A/B 


1.30 


ns 


1.37 


ns 


-1.14 


ns 


-1.20 


ns 



fad: Fatty acyl desaturase (A5 and A6 activities); NADH1: NADH dehydrogenase subunit 1; UCP2: Uncoupling protein 2; elovl: MTHFD1: Methylenetetrahydrofolate 
dehydrogenase (NADP+ dependent) 1 -like; PA2G4: Proliferation-associated 2G4b; RP60S: 60S acidic ribosomal protein; MYO: Myosin heavy chain; PCNA: 
Proliferating cell nuclear antigen; CYP1A: Cytochrome P4501A; ABCA1: ATP-binding cassette sub-family A member 1; PSMB8: Proteasome subunit beta type-8 
precursor; TAGLN: Transgelin 2; COL1A2: Collagen type I alpha 2; elovl2: Elongase 2; FAS: Fatty acid synthase; ACO: Acyl-CoA oxidase; CPT1: Carnitine 
palmitoyltransferase I; PPARs; Peroxisome prol iterator-activated receptors; GST: glutathione S-transferase; MT-A: Metallothionein A; CAT: catalase; SOD: superoxide 
dismutase; CASP; caspase. 
*Only marginally non-significant. 

Genes correspond to features significantly affected by either diet (D) or genotype (G) in the microarray analysis or involved in biological processes being analysed 
in further detail. Values are normalized (by cofilin-2 and elf-1 a) gene expression ratios between fish fed VO in relation to FO for each family group or of Lean fish 
in relation to the Fat group when fed either one of the diets. Values in bold are significantly different, at p<0.05 (REST 2008); ns - non-significant. 
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collagen (COL1A2), while transgelin 2 (TAGLN) expres- 
sion did not differ between family groups, and for ATP- 
binding cassette sub-family A member 1 (ABCA1) there 
was an inverse change in expression. 

As well as validation above, RT-qPCR was used to fur- 
ther analyse genes/pathways identified by microarray and 
published data as potentially interesting, including lipid 
metabolism, xenobiotic and oxidative stress, and apoptosis. 
One was LC-PUFA biosynthesis, given that A5fad was sig- 
nificantly affected by diet in the microarray analysis, with a 
stronger response in Fat fish, whereas A6fad showed a sig- 
nificant diet x genotype interaction (supplementary file 1) 
confirmed by RT-qPCR. The A6fad transcript was only sig- 
nificantly up-regulated in Fat fish fed VO, compared to 
FO, and Lean fish showed higher levels of A6fad expres- 
sion than Fat fish when fed FO, while the opposite trend 
was noted when fish were fed VO. Fatty acyl elongases 
(elovls) were also quantified and their expression (except 
elovl5a) broadly followed that of fads: significantly up- 
regulated when dietary VO replaced FO in the Fat group. 
Additionally, elovl5b and, particularly, elovl2 showed a 
trend for increased expression (high fold-change) in Lean 
fish, compared to Fat fish, when fed FO, while an opposite 



trend was observed when salmon were fed VO. Although 
genes involved in fatty acid synthesis and oxidation showed 
few significant differences, expression of fatty acid synthase 
(FAS) was up-regulated in fish fed VO, but only in Lean 
fish. The expression of peroxisome proliferator-activated 
receptors (PPARs), involved in the regulation of multiple 
lipid metabolism genes, was determined but only PPARy 
showed any significant change, being up-regulated in the 
Fat group when dietary VO replaced FO. Of the xenobiotic 
and oxidative metabolism genes assayed, apart from CYP1A, 
only catalase (CAT) was affected by diet (up-regulated in 
fish fed VO) and only significantly in the Lean family. In 
contrast, metallothionein A (MT-A) showed higher expres- 
sion in Fat fish, but only when fed VO, while a marginal 
down-regulation was observed when comparing VO and 
FO-fed fish in the Lean group. Of genes related to apoptosis, 
CASP3B was up-regulated by VO in Lean fish whereas a 
similar fold-change was marginally non-significant in the 
Fat fish. 

Intestine fatty acid composition 

The levels of most fatty acids in pyloric caeca were 
affected by diet, whereas genotype had no significant 



Table 4 Fatty acid compositions (% of total fatty acids) and lipid contents (% of wet weight) of intestinal tissue of 
Lean and Fat family groups fed diets containing either FO or VO 



FO VO Two-way ANOVA 





Fat 


Lean 


Fat 


Lean 


D 


G 


DxG 


Fatty acid 


















Total saturated 


23.70 ± 


1.21 


22.82 ± 1 .09 


18.57 ± 0.73 


20.08 ± 1 .38 


0.000 


0.742 


0.066 


Total monoenes 


38.25 ± 


7.15 


44.45 ±5.31 


44.98 ± 2.38 


39.34 ± 4.49 


0.717 


0.928 


0.046 


18:2n-6 


4.19 ± 


0.67 


4.51 ± 0.24 


13.14 ± 0.59 


11.28 ± 1.12 


0.000 


0.139 


0.041 


18:3n-6 


0.10 ± 


0.02 


0.10 ± 0.02 


0.28 ± 0.12 


0.25 ± 0.13 


0.003 


0.611 


0.611 


20:2n-6 


0.57 ± 


0.06 


0.58 ± 0.07 


1 .54 ± 0.40 


1 .65 ± 0.06 


0.000 


0.199 


0.510 


20:3n-6 


0.17 ± 


0.06 


0.14 ± 0.02 


1.01 ± 0.12 


1 .47 ± 0.26 


0.000 


0.056 


0.019 


20:4n-6 


1.20 ± 


0.60 


0.85 ± 0.31 


0.74 + 0.18 


1 .40 ± 0.29 


0.785 


0.419 


0.016 


22:4n-6 


0.08 ± 


0.06 


0.05 ± 0.06 


0.00 ± 0.00 


0.09 ± 0.07 


0.337 


0.337 


0.038 


22:5n-6 


0.21 ± 


0.06 


0.16 ± 0.07 


0.07 ± 0.02 


0.12 ± 0.05 


0.002 


1.000 


0.128 


Total n-6 PUFA 


6.51 ± 


0.18 


6.38 ± 0.36 


16.79 ± 0.86 


16.27 ± 0.67 


0.000 


0.410 


0.781 


18:3n-3 


1.46 ± 


0.38 


1.67 ± 0.21 


4.79 ± 0.53 


3.71 ± 0.75 


0.000 


0.211 


0.021 


18:4n-3 


1.36 ± 


0.60 


1 .67 ± 0.66 


0.85 ± 0.31 


0.55 ± 0.20 


0.002 


0.712 


0.127 


20:3n-3 


0.27 ± 


0.05 


0.33 ± 0.14 


0.67 ± 0.18 


0.62 ± 0.04 


0.000 


0.768 


0.384 


20:4n-3 


1.77 ± 


0.28 


2.16 ± 0.46 


0.97 ± 0.13 


0.83 ± 0.10 


0.000 


0.504 


0.061 


20:5n-3 


6.55 ± 


0.55 


5.79 ± 0.62 


3.60 ± 0.90 


4.40 ± 0.72 


0.000 


0.682 


0.058 


22:5n-3 


2.82 ± 


0.81 


2.30 ± 0.74 


1.57 ± 0.36 


2.54 ± 0.77 


0.157 


0.530 


0.052 


22:6n-3 


16.68 ± 


6.68 


11.65 ± 5.65 


6.97 ± 1.58 


11.50 ±3.10 


0.057 


0.943 


0.050 


Total n-3 PUFA 


30.90 ± 


5.94 


25.56 ± 5.19 


19.41 ± 2.18 


24.15 ± 3.76 


0.015 


0.962 


0.042 


Total PUFA 


38.06 ± 


5.96 


32.73 ± 5.42 


36.45 ± 2.08 


40.59 ± 3.16 


0.206 


0.794 


0.059 


Total lipids 


7.50 ± 


3.75 


9.56 ± 3.81 


6.97 ± 1 .07 


4.59 ± 1 .05 


0.081 


0.804 


0.111 



Results are means ± SD (n = 4) and p-values of two-way ANOVA are presented for diet (D), genotype (G) and interaction between both factors. 
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effect (Table 4). However, some fatty acids also showed a 
significant diet x genotype interaction, indicating that 
the effect of diet depended on the genetic background of 
the fish. For instance, interactions were observed for 
some LC-PUFA as a result of higher levels being found 
in the Lean group, compared to Fat, when fish were 
fed VO, while the reverse was observed when fed FO. 
Another unexpected result was that similar levels of 
DHA in FO- and VO-fed Lean fish meant that, in 
spite of substantial differences in Fat fish fed the two 
diets, the effect of diet on DHA was marginally non- 
significant. Similarly, levels of EPA and 22:5n-3 be- 
tween FO- and VO-fed fish were noticeably closer in 
the Lean group. 



Proteomic analysis 

Of the protein spots identified as being differentially 
expressed between diets or genotypes (Additional files 2 
and 3), only 17 and 29 could be excised and, of these, 
only 9 and 20, respectively, returned reliable identifica- 
tions by peptide fragment fingerprinting (Figure 1). Pro- 
teins significantly up-regulated by dietary VO are likely 
implicated in xenobiotic/drug metabolism (epoxide 
hydrolase 2, EPHX2), protection from oxidative stress 
(hemopexin-like protein, HPX and peroxiredoxin-1, 
PRDX1) and induction of apoptosis and inflammatory 
responses (galectin-2, LGALS2). Those proteins down- 
regulated by dietary VO included proteins responsible for 
protein folding and involved in signalling (2-peptidylprolyl 



CI 




Figure 1 Proteins significantly affected by diet or genotype. Cropped image of a scanned 2-D polyacrylamide gel showing all spots 
identified as being differentially expressed (P<0.05, two-way ANOVA) in fish fed either FO or VO (underlined spot numbers), or between family 
groups (remaining spot numbers); the three underlined spot numbers containing an * were differentially affected by both diet and genotype. 
Spots that returned reliable identifications by peptide fragment fingerprinting are identified by their abbreviation. 



Table 5 Proteins differentially regulated by diet 



Spot No 


Protein ID 


Accession No. 
(NCBInr) 


Theoretical Mw 
(KDa)/pl 


Protein 
Score 3 


No. Matched peptides 
(MS/MS) b 


Best score peptide 


VO/FO Lean 


VO/FO Fat 


p-value 


1136 


Epoxide hydrolase 2 (EPHX2) 


gb|ACI331 29.1 1 


54.2/5.52 


84 


1 


GGLFVGLPDEIPR 


1.23 


1.40 


0.0190 


1148 


Hemopexin-like protein* (HPX) 


emb|CAA92147.1| 


50.4/5.61 


169 


1 


VHLDAITSDDAGNIYAFR 


1.32 


1.56 


0.0002 


1151 


Hemopexin-like protein* (HPX) 


ref|NP_001 10461 7.1 1 


51.0/6.18 


161 


1 


VHLDAITSDDAGNIYAFR 


1.38 


1.54 


0.0001 


2683 


Peroxiredoxin 1 (PRDX1) 


gb|ACI67145.1| 


22.0/6.42 


757 


8 


SISTDYGVLKEDEGIAYR 


1.14 


1.09 


0.0310 


3168 


2-peptidylprolyl isomerase A 
(PPIA) 


ref|NP_001 135161.1 1 


17.5/7.61 


101 


1 


VYFDITIGDTPAGR 


-1.46 


-1.18 


0.0170 


3243 


Myosin light chain smooth 
muscle isoform (MYL) 


ref|NP_998803.1| 


16.9/4.47 


100 


2 


EAFLLFDR 


-1.23 


-1.11 


0.0490 


3331 


similar to H2A histone family, 
member J (H2A) 


ref|XP_001521566.1| 


13.3/10.84 


74 


1 


AGLQFPVGR 


-1.58 


-1.18 


0.0250 


3334 


Histone H2A** (H2A) 


emb|CAA25528.1| 


13.7/10.88 


37 


1 


AGLQFPVGR 


-1.48 


-1.18 


0.0210 


3445 


Galectin 2 (LGALS2) 


gb|ACN10131.1| 


14.8/5.93 


314 


2 


SGASSFSINIGHDSDNYALHFNPR 


1.17 


1.15 


0.0170 



a The protein score probability limit (where P < 0.05) is 73. 

b Peptides with confidence interval above 95% were considered. 

* Equivalent to warm-temperature-acclimation-related-65 kDa-protein. 

**Same protein identification obtained by MS/PMF (peptide mass fingerprinting): Score = 76; 7 peptides matched. 

Only reliable identifications of Actinopterygii, obtained by mass spectrometry (MALDI-TOF-MS/MS) analysis and searches in MASCOT, NCBI and ExPASy (Mw/pl) databases are shown. Expression ratios between fish fed 
VO and FO diets are given for each family group, as well as p-value for diet (two-way ANOVA, DeCyder V7.0). 



Table 6 Proteins differentially regulated by genotype 



Spot No Protein ID Accession No. Theoretical Mw Protein No. Matched peptides Best score peptide Lean/Fat FO Lean/Fat VO p-value 







(NCBInr) 


(KDa)/pl 


Score 3 


(MS/MS) b 










1148 


Hemopexin-like protein* (HPX) 


emb|CAA92147.1| 


50.4/5.61 


169 


1 


VH L D A ITS D D AG N I Y A F R 


-1.09 


-1.29 


0.0270 


1389 


Calreticulin precursor (CALR) 


gb|ACI32936.1| 


47.6/4.33 


81 


2 


FEPFSNEGK 


-1.28 


-1.33 


0.0026 


1441 


Alpha-enolase (EN01) 


ref|NP_001 133366.1 1 


47.0/5.91 


405 


3 


AAVPSGASTGIYEALELR 


1.28 


1.2 


0.0066 


1714 


Heat shock protein 70 (HSP70) 


ref|NP_990334.1| 


70.8/5.47 


1 1 1 


1 


IINEPTAAAIAYGLDKK 


1.16 


1.77 


0.0260 


2134 


Heterogeneous nuclear 
ribonucleoprotein AO (HNRNPAO) 


gb|ACI67551.1| 


28.8/9.10 


235 


2 


LFVGGLNVDTDDDGLRK 


-1.19 


-1.09 


0.0460 


2154 


Calponin 1 (CNN1) 


ref|NP_001 139857.1 1 


33.2/8.56 


565 


4 


KINTSPQNWHQLENIGNFVR 


-1.34 


-1.51 


0.0350 


2157 


Calponin 1*** (CNN1) 


ref|NP_001 139857.1 1 


33.2/8.56 


43 


1 


YDPQKEEELR 


-1.34 


-1.47 


0.0240 


2213 


Caspase 3 (CASP3) 


ref|NP_001 133393.1 1 


31.0/5.97 


254 


2 


I PVEADFLYAYSTAPGYYSWR 


-1.15 


-1.24 


0.0300 


2215 


Caspase 3 (CASP3) 


ref|NP_001 133393.1 1 


31.0/5.97 


279 


2 


VANDQWQQIQQLLSK 


-1.14 


-1.31 


0.0130 


2282 


Voltage-dependent anion 
channel 2-2 (VDAC2) 


gb|ACH71030.1| 


30.1/8.85 


214 


1 


VN N N S LVGVG^TQTLR P G VK 


-1.24 


-1.19 


0.0160 


2458 


Annexin A4 (ANXA4) 


gb|ACI69495.1| 


28.4/5.22 


218 


3 


NHLLQVFK 


-1.14 


-1.11 


0.0300 


2487 


Dihydropteridine reductase 
(DHPR) 


gb|ACI67281.1| 


15.7/8.46 


81 


1 


QSWWSTISSHLATR 


-1.14 


-1.11 


0.0300 


2499 


Pancreatic alpha-amylase 
precursor (AMY2) 


ref[NP_001036176.1| 


57.4/6.89 


86 


1 


ALVFVDNHDNQR 


-1.37 


-1.27 


0.0021 


2522 


Triosephosphate isomerase 
1 b (TP 1 1 ) 


ref|NP_001 1331 74.1 1 


26.6/7.63 


336 


2 


LDPNTEWCGAPSIYLEFAR 


1.26 


1.16 


0.0009 


2527 


Histone cluster 1 (H2A) 


ref|NP_001 086775.1 1 


14.0/10.88 


94 


1 


AGLQFPVGR 


1.12 


1.11 


0.0074 


2677 


Transgelin (TAGLN) 


gb|ACM09025.1| 


21.7/7.69 


219 


1 


DGCVLSELINSLHK 


-1.61 


-1.35 


0.0350 


2717 


Proteasome beta 1 subunit 
(PSMB1) 


ref|NP_001 003889.1 1 


26.1/6.32 


87 


1 


GAVYSFDPVGSYQR 


-1.29 


-1.12 


0.0140 


2903 


Peroxiredoxin-1(PRDX1) 


gb[ACI67145.1| 


22.0/6.42 


115 


1 


QITINDLPVGR 


-1.62 


-1.99 


0.0069 


3345 


Retinol-binding protein II, 
cellular (RBP2) 


ref|NP_001 139954.1 1 


15.6/5.44 


70 


2 


AIDIDFATR 


-1.32 


-1.1 


0.0330 


3457 


Alpha globin (HBA) 


emb|CAA65949.1| 


15.1/9.19 


170 


2 


TYFSHWADLSPGSAPVK 


1.03 


2.09 


0.0320 



a The protein score probability limit (where P < 0.05) is 73. 

b Peptides with confidence interval above 95% were considered. 

* Equivalent to warm-temperature-acclimation-related-65 kDa-protein. 

*** Same protein identification obtained by MS/PMF (peptide mass fingerprinting): Score = 193; 27 peptides matched. 

Only reliable identifications of Actinopterygii, obtained by mass spectrometry (MALDI-TOF-MS/MS) analysis and searches in MASCOT, NCBI and ExPASy (Mw/pl) databases are shown. Expression ratios between Lean 
and Fat family groups are given for each dietary treatment, as well as p-value for genotype (two-way ANOVA, DeCyder V7.0). 
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isomerase A, PPIA), actin-based motility (myosin light 
chain, MYL) and DNA replication, repair or transcription 
(histone H2A, H2A) (Table 5). Proteins affected by geno- 
type encompass a variety of pathways, of which only a few 
are related to metabolism, namely carbohydrate (alpha- 
enolase, ENOl; pancreatic alpha-amylase precursor, 
AMY2; triosephosphate isomerase lb, TPI1), folate (dihy- 
dropteridine reductase, DHPR) or retinol (retinol-binding 
protein II, RBP2) metabolism (Table 6). Other proteins 
may have potentially multiple roles but can broadly be 
assigned roles in response to oxidative and cellular stress 
(HPX; PRDX; heat shock protein 70, HSP70), oxygen 
transport (alpha globin, HBA), signal transduction (calre- 
ticulin, CALR), transcription/RNA repair (heterogeneous 
nuclear ribonucleoprotein AO, HNRNPAO; histone H2A, 
H2A), apoptosis (CASP3), cellular transport, potentially 
also associated with apoptosis (voltage-dependent anion 
channel 2-2, VDAC2; annexin A4, ANX4), and proteolysis 
(proteasome beta 1 subunit, PSMB1). As with the micro- 
array analysis, a few proteins with a more structural func- 
tion and particularly associated with tissue contractile 
properties were affected by genotype, showing lower levels 
in Lean fish. These included calponin-1 (CNN1) and 
transgelin (TAGLN), the latter which was also found to be 
significantly affected by microarray, albeit up-regulated in 
Lean fish. Most proteins significantly affected by genotype 
showed lower levels of expression in the Lean group, with 
the exception of ENOl, HSP70, TPI1, H2A and HBA. 

Discussion 

Dietary plant ingredients can induce chronic intestinal 
inflammatory conditions in salmonids that can ultim- 
ately result in carcinogenesis [25]. This extreme reaction 
is rare and usually associated with soy protein at high 
levels [26]. Dietary n-3 LC-PUFA have important anti- 
inflammatory and anti-carcinogenic effects in mamma- 
lian intestine [23]. Therefore, use of feeds containing 
high percentages of plant proteins combined with re- 
placement of FO by VO, as is now prevalent in the in- 
dustry, requires studies on dietary effects on intestinal 
transcriptomes and proteomes. However, interpretation 
of the data was difficult as the effects on dietary treat- 
ments and/or family groups were subtle, as also 
observed in liver transcriptome [9], and is typical of this 
type of experiment [6,8,11,27]. Partly as a consequence, 
validation of the microarray data gave variable results, 
from perfect match, to opposite changes in a few, al- 
though effects observed in the microarray, with fold 
changes as low as 1.2 (e.g., NADH1 and PA2G4) were 
validated by RT-qPCR. In view of the whole genome du- 
plication event that occurred in salmonids [28], gene ex- 
pression studies are often more challenging due to the 
presence of highly similar genes which may hybridize 
with cDNA probes presenting low specificity, further 



complicated if similar transcripts, corresponding to 
duplicated genes, are differentially regulated [10]. None- 
theless, the presence of several features related to spe- 
cific processes in both the transcriptomic and proteomic 
analysis gave supporting evidence to the pathways likely 
differentially affected by dietary oil and genetic back- 
ground related to flesh adiposity. 

Intestinal LC-PUFA biosynthesis capacity is differentially 
affected by diet and genotype 

Considering whether genetic selection for fish families 
showing better adaptation to more sustainable feeds 
might be a viable approach to develop aquaculture, one 
outcome of this investigation was to establish if effects 
of diet on expression of LC-PUFA biosynthesis genes 
depended on genotype, as shown in the liver transcrip- 
tome of these fish [9]. This was not seen in the hepatic 
transcriptome of European sea bass families showing dif- 
ferent growth rates when fed a vegetable diet but, in this, 
case similar LC-PUFA profiles were also noted in both 
genotypes in response to the vegetable diet [11]. In both 
salmon tissues, differences in n-3 LC-PUFA content be- 
tween fish fed FO or VO were smaller in the Lean family 
group. This was due to higher levels of n-3 LC-PUFA in 
Lean salmon, compared to Fat, when fish were fed VO, 
but higher amounts in the Fat family group when fed 
FO. However in liver, up-regulation of LC-PUFA biosyn- 
thesis when fish were fed VO was much larger in the 
Lean family group, whereas in intestine the same indivi- 
duals only showed significant up-regulation in Fat fish. 
This appears contradictory but can be explained by the 
differential tissue n-3 LC-PUFA contents. Although the 
difference was smaller in Lean fish compared to Fat fish 
in both tissues, in liver there was still a considerable dif- 
ference in n-3 LC-PUFA levels between fish fed FO or 
VO, while in intestine levels were similar. PUFA have 
important activities on transcription factors, either as 
direct ligands or through effects on membrane compos- 
ition [29], affecting transcription of many genes involved 
in lipid metabolism, including desaturases and elongases 
[30,31]. In salmon, regulation of genes of LC-PUFA bio- 
synthesis that are known to respond to dietary compos- 
ition, i.e., A5fad, A6fad, elovl5b and elovl2 [32-34], 
appear to show high plasticity and are likely under feed- 
back regulation by tissue n-3 LC-PUFA. Both studies 
suggest that the Lean family group may show an 
enhanced response to low dietary n-3 LC-PUFA, with 
greater up-regulation of biosynthesis when fed VO. In 
contrast to liver, this response was sufficient in intestine 
to maintain tissue n-3 LC-PUFA, particularly DHA, at 
similar levels to fish fed FO. Considering that differences 
in desaturase expression between the Fat and Lean fish 
were only significant when FO but not VO, was fed, sug- 
gests that the likely mechanism is through negative 
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feedback by high levels of n-3 LC-PUFA rather than 
positive feedback from low levels of LC-PUFA and/or 
higher levels of shorter chain precursors. 

Other dietary effects on lipid metabolism 

Transcriptional regulation of desaturases and elongases 
by LC-PUFA may involve both PPARa and sterol regula- 
tory element binding protein lc (SREBP-lc) [30,31]- In 
liver, expression of A5fad, A6fad, elovl2, PPARa, and 
possibly PPARp, appeared co-ordinately regulated by 
diet depending on genotype, while PPARy was not 
affected [9]. In intestine, however, expression of PPARa 
and PPARp was not affected by either diet or genotype, 
while PPARy was up-regulated by dietary VO, signifi- 
cantly in Fat fish. This suggests that dietary regulation of 
lipid metabolism genes in fish intestine might differ to 
mammals, where PPARa showed differential expression 
in response to dietary EPA and DHA in murine intestine 
[35]. Reasons for differential regulation of PPARs be- 
tween salmon liver and intestine are unclear, but may be 
due to different patterns of tissue expression. In plaice 
and seabream, there was no nutritional regulation (fast- 
ing/feeding) of PPARs in the intestine, where PPARy 
was the dominant isotype, in contrast to liver where 
PPARa was dominant [36]. PPARy in both mammals 
and fish is predominantly expressed in adipose tissue 
and promotes adipocyte differentiation and lipid storage 
[36-38]. In mammals, PPARy activates the expression of 
genes characteristic of mature adipocytes and adipogen- 
esis, including FAS [37,38] and hence the expression of 
PPARy, up-regulated in salmon fed VO, might be related 
to increased expression of FAS. However, increased 
PPARy expression was only significant in Fat fish 
whereas FAS was significantly up-regulated only in Lean 
salmon. As fish PPARy is functionally the most different 
of the three isotypes compared to mammalian PPARs, 
and is expressed more widely in fish tissues that in 
mammals, other mechanisms and functions may under- 
lie the observed regulation [35]. 

In this study, the hypotriglyceridemic effect of LC- 
PUFA, well established in mammals [39,40], was also 
observed in salmon intestine. Lipogenesis was down- 
regulated in FO-fed fish, as demonstrated by decreased 
FAS expression (RT-qPCR) and the presence of a tran- 
script containing a beta-ketoacyl synthase domain, a 
component of FAS (microarray). The differences in FAS 
expression were not as marked as in liver [9] and were 
only significant in Lean fish but, together with the LC- 
PUFA biosynthesis data, demonstrate the active role of 
salmon intestine in lipid metabolism [21]. However, des- 
pite up-regulation of lipogenesis by dietary VO, lipid ac- 
cumulation in enterocytes was lower than in fish fed FO, 
contrary to previous reports of VO promoting lipid ac- 
cumulation in enterocytes [16-18]. 



In contrast, the hypotriglyceridemic effect of LC-PUFA 
did not involve the typical increase in p-oxidation [39], 
reported in mice intestine [36,40] . As in liver, no changes 
were observed in the expression of p-oxidation genes car- 
nitine palmitoyltransferase I (CPT1) and acyl-CoA oxi- 
dase (ACO) [10]. Nonetheless, effects of dietary lipid on 
energy metabolism were observed in intestine. In particu- 
lar, UCP and transcripts involved in the mitochondrial 
electron transport chain, including components of cyto- 
chrome c oxidase, NADH1 and ubiquinol- cytochrome c 
reductase complexes, and the mitochondrial metabolite 
transporter SCaMC-2, were slightly down-regulated by 
dietary VO, possibly suggesting reduced energetic metab- 
olism. EPA may act as a mitochondrial proliferator in 
both rat and salmon liver [41,42], which might also ex- 
plain this result. 

Vegetable oils as potential sources of contaminants in 
alternative diets for aquaculture? 

The use of FO in feeds has also been questioned in rela- 
tion to levels of persistent organic pollutants, POPs [43], 
whereas VO are generally considered safer alternatives to 
dietary FO [44]. However, results suggested that VO had 
greater impact on xenobiotic metabolism in salmon in- 
testine than FO, with the transcriptome and proteome 
both showing up-regulation of transcripts and proteins 
involved in detoxification and protection from oxidative 
stress in fish fed VO. This was surprising, considering ex- 
pression of these genes has been associated with FO sup- 
plementation and higher levels of organic contaminants 
and/or increased peroxidative susceptibility of LC-PUFA 
[22,45,46]. In particular, and linked to a detoxification 
role, up-regulation of CYP1A transcript and epoxide 
hydrolase 2 (EPHX2) protein was found in fish fed VO. 
Cytochrome P450 1A metabolizes many exogenous and 
endogenous molecules, including pollutants (e.g., poly- 
chlorinated biphenyls-PCBs, polycyclic aromatic hydro- 
carbons -PAHs, dioxins and dibenzofurans) as well as 
several metabolic products (e.g., steroids, bile acids). Its 
expression is sensitive to low levels of contaminants and 
therefore it is a commonly used marker [47,48]. Lipid 
peroxidation and antioxidant enzymes are also biomar- 
kers for environmental xenobiotic contamination in fish, 
as the catalytic actions of detoxifying enzymes like 
CYP1A produce high levels of reactive oxygen species 
(ROS) [49,50]. We observed up-regulation of CAT and of 
a selenoprotein transcript, as well as of HPX and PRDX1 
proteins in salmon fed VO. CAT and PRDX1 catalyze the 
decomposition of hydrogen peroxide [51], hemopexin 
prevents heme-mediated oxidative stress [52], and sele- 
noproteins have diverse roles, including selenium trans- 
port and antioxidant defense [53]. Expression of other 
genes involved in protection from oxidative stress like 
GST and SOD were not affected, but MT-A was down- 
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regulated in the Lean family group fed VO. MT-A, be- 
sides having an antioxidant role, also protects against 
heavy metal toxicity and maintains physiological zinc 
homeostasis [54] and hence could be responding to con- 
taminants more abundant in FO (e.g., heavy metals) [55]. 

Previous analysis of contaminants in comparable feeds 
using the same standard FO and a similar VO blend, 
showed that levels of several POPs, including organo- 
chlorine pesticides, and heavy metals were substantially 
lower, but polycyclic aromatic hydrocarbons (PAH) were 
10-fold increased in the VO diet compared to the FO 
diet [55]. PAH derive from incomplete combustion of 
organic compounds and can be found in high concentra- 
tions in fats, including VO, through multiple routes of 
contamination, before, during or after oil processing 
[56]. They are metabolized by both CYP1A1 and 
EPHX2, among other enzymes, and CYP1A1 is induced 
by PAH in mammals [57]. Therefore, higher PAH levels 
in the VO diet might explain, at least partly, the results 
obtained in the present study although, unlike POPs, 
PAHs are not persistent and are readily eliminated from 
fish tissues [58]. 

High doses of PAH result in significant intestinal 
hyperplasia in fish, with an increase in cell proliferation 
and faster epithelial turnover [59]. Previous studies on 
intestinal gene expression in fish indicated a reduction 
in cell proliferation or differentiation associated with 
dietary FO replacement by VO [21,22], possibly due to 
lower levels of membrane LC-PUFA and reduced oxida- 
tive stress [22]. In the present study, no major impact on 
cell proliferation was apparent in the intestinal transcrip- 
tome or proteome data. Two transcripts related to cell 
proliferation, PA2G4 and cyclin Gl, were slightly down- 
regulated in fish fed VO, but in mammals these have op- 
posing effects [60,61] and, furthermore, two mammalian 
PA2G4 isoforms have been shown to have opposite 
effects in cellular proliferation [62] and hence results are 
inconclusive. 

Previously, expression of caspases, effectors of con- 
trolled cell death or apoptosis, was affected by replace- 
ment of dietary FO by VO in fish [21,22]. Apoptosis is 
particularly important in organs with high rates of cellu- 
lar turnover such as intestine but, in addition to main- 
taining normal gut function, apoptosis may be affected by 
pathological or toxic conditions, including those induced 
by environmental chemical contaminants [63,64], In the 
present study, expression of CASP3B was up-regulated 
in salmon fed VO, particularly in the Lean family group 
and a similar, non-significant trend was observed for 
CASP6A/B. As ROS are important signalling molecules 
in apoptotic processes, these results could be linked to a 
cytotoxic effect causing increased oxidative stress in VO 
[64]. Relevant to the above was the up-regulation of 
galectin 2 (LGALS2) in the proteome of salmon fed VO. 



Galectins are pleiotropic regulators of immune functions 
and are up-regulated by injury and infectious conditions, 
have well-recognized modulatory roles in mammalian 
intestinal inflammatory diseases, and their mode of ac- 
tion involves induction of apoptosis [65,66]. The lack of 
major effects on cell proliferation and only slight up- 
regulation of CASP3 and LGALS2 suggests that any 
contaminant doses experienced by the fish were unlikely 
to have caused any serious morphophysiological damage 
in the intestine. As similar trends were not seen in the 
hepatic transcriptome of these individuals [9], this may 
suggest intestine can potentially metabolize and detoxify 
xenobiotics present in the diet. Furthermore, there were 
no growth or general performance issues with these fish 
[5]. Therefore, the data do not imply abnormal gastro- 
intestinal functions or effects on final product quality. 

Effect of genotype in intestinal transcriptome 
and proteome 

Contrary to diet, genotype did not have a major impact 
on metabolism genes, apart from transcripts related to 
the proteasomal degradation pathway including a strong 
down-regulation of PSMB8 in Lean fish, particularly fed 
VO. This gene has been recently found to have a mo- 
lecular evolution history that suggests a very strong se- 
lective pressure for its functional dimorphism to be 
maintained in vertebrates [67]. Two different alleles, A- 
type and F-type, can be found in basal vertebrate species, 
including Atlantic salmon. The PSMB8F lineage was lost 
in common ancestors of higher teleosts and tetrapods 
but was then independently revived de novo through the 
appearance of F-type alleles within the PSMB8A lineage. 
In this study we did not find evidence of significant dif- 
ferences between families groups for the A-type allele in 
the transcriptomic analysis as the probe showing signifi- 
cant variation between families in the microarray corre- 
sponded to the PSMB8F allele. Hence, it was also the 
F-type transcript that was validated by RT-qPCR, using 
type-F specific primers. However, further studies would 
be required to confirm this and to assess the functional 
significance of this result. On the other hand, expression 
of PSMB1 was down-regulated in the intestine proteome 
of Lean fish. Proteolysis through this pathway is essential 
for many cellular processes, including the cell cycle, sig- 
nalling, cellular defence and responses to oxidative stress 
[68]. Therefore, this response might be related to de- 
fence against cellular stress, as another difference be- 
tween the two family groups was related to xenobiotic 
and oxidant metabolism. Apart from lower expression of 
a CYP1A transcript in Lean fish, two proteins with anti- 
oxidant roles, HPX and PRDX, were down-regulated in 
the proteome of Lean compared to Fat fish. Alpha glo- 
bin, or haemoglobin alpha (HBA), a major component of 
blood and potent mediator of oxidative stress, can have 
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both protective and damaging effects depending on 
complex interactions in H 2 0 2 -rich environments [69]. 
However, given its opposite regulation to HPX, whose 
main role is to scavenge heme and protect from its toxic 
effects [52], up-regulation of HBA in Lean fish may indi- 
cate heme-mediated oxidative stress. 

The apoptotic pathway may be differentially affected 
by genotype, with down-regulation of CASP3, VDAC2 
and ANXA4 in the Lean family group, the latter two 
transport proteins having well-recognized roles in apop- 
tosis [70,71]. In contrast, heat shock proteins that pro- 
tect against environmental stresses were increased in the 
intestine transcriptome and proteome of Lean salmon. 
This response could be associated with the observed 
changes in the ubiquitin-proteasome degradation sys- 
tem, as the systems have been functionally coupled in 
mammals. Thus, moderate exposure to a heat shock can 
cause a transient increase in intracellular proteolysis by 
the ubiquitin-proteasome pathway, followed by a phase 
of slower or even inhibited protein degradation [72,73]. 
Furthermore, Pirkkala et al. [74] demonstrated transcrip- 
tional induction of heat shock genes when proteasome 
activity was down-regulated. However, judging by the 
fold-changes, these effects are only relevant when fish 
were fed VO, and hence could be more related to dietary 
changes. Collectively, the data may indicate higher sensi- 
tivity of Lean fish to environmental or endogenous stres- 
ses due to replacement of dietary FO by VO. 

The predominant influence of genotype was in the ex- 
pression of intestinal transcripts of structural proteins, 
particularly collagen alpha chains, but also osteonectin, 
TAGLN, troponin I and keratocan, which were up- 
regulated in Lean fish, whereas troponin C was down- 
regulated. Furthermore, CNN1 and TAGLN were 
down-regulated in the intestinal proteome in Lean fish. 
Collagen, the main component of connective tissue, 
helps to maintain the structural integrity of tissues [75], 
while osteonectin is an extracellular matrix glycoprotein 
with high affinity towards collagen and whose expression 
has been associated with remodelling processes in tis- 
sues, including human intestine during development/ 
morphogenesis and in diseased mucosa [76]. Troponin, 
TAGLN and CNN1 are all involved in actin binding, 
actin-myosin interaction and muscle contraction. The 
inverse regulation of troponins is not conflicting as they 
have different roles in actomyosin cross-bridge forma- 
tion and contraction; binding of troponin C to Ca 2+ 
induces conformational changes that counteract the in- 
hibitory action of troponin I [77]. Expression of TAGLN 
transcript and protein showed opposite effects but a lack 
of correlation between transcriptomic and proteomic 
data is not unprecedented [78]. As discussed above, this 
result might also be explained by the presence of similar 
duplicated genes in Atlantic salmon that are regulated 



differently. Transcriptomic results were validated by RT- 
qPCR for COL1A2, although only significantly when fish 
were fed the VO diet, for which fold changes were 
higher. In addition, in the microarray results differences 
in expression of structural proteins between family 
groups were consistently more accentuated in fish fed 
VO which could suggest a cumulative effect of diet. Fur- 
thermore, MYO was up-regulated in fish fed VO com- 
pared to FO but only in Lean fish (RT-qPCR), and 
significant diet x genotype interactions were found for 
alpha-actinin 1, tubulin beta-2 chain and procollagen- 
lysine 2-oxoglutarate 5-dioxygenase 2 (Additional file 1), 
which were up-regulated in Lean salmon, compared to 
Fat, but only when fed VO. In cod, replacement of FO 
by VO resulted in changes in intestinal expression of 
structural genes with the potential to alter the structural 
and mechanical properties of the intestinal muscle layer, 
including a range of actin-binding transcripts [21]. 

The present study is the first investigation of the influ- 
ence of genetic background of families with different 
flesh adiposity phenotypes on intestinal gene expression 
of a fish species. Effects were subtle and consequently 
their potential impacts difficult to fully assess. However, 
if genetic selection for families better adapted to alterna- 
tive formulations is an approach taken in the future, the 
potential for genotype-specific differences being exacer- 
bated when VO replaces dietary FO should be further 
examined to assess the consequences of these changes in 
intestinal gene expression. 

Conclusions 

Metabolic activity, particularly lipid and energy, of intes- 
tinal tissue responded to dietary lipid composition but 
was also affected by genotype. The LC-PUFA biosyn- 
thesis pathway, typically up-regulated when salmon are 
fed VO, was especially influenced by genetic back- 
ground. The Lean fish showed an enhanced response to 
low dietary n-3 LC-PUFA and the expression of A5fad, 
A6fad, elovl5b and elovl2 in the intestine showed high 
plasticity and was reflected in tissue biochemical com- 
position indicating that their transcriptional regulation 
might be under feedback control by n-3 LC-PUFA, 
mainly DHA. Lower n-3 LC-PUFA in VO increased lipo- 
genesis in Lean salmon, assessed by expression of FAS, 
while p-oxidation appeared unaffected, although tran- 
scripts involved in mitochondrial respiratory or electron 
transport chains were down-regulated, suggesting 
reduced activity in fish fed VO. Higher expression of 
genes and proteins involved in xenobiotic metabolism 
(CYP1A and EPHX2), antioxidant defence (CAT, HPX 
and PRDX1), and apoptosis (Casp3B) were observed in 
VO-fed fish, suggesting they might be responding to 
higher levels of contaminants, particularly PAH, in the 
diet. However, the intestine appeared able to metabolize 
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and detoxify xenobiotic substances potentially present in 
the diet without major deleterious effects. Nonetheless, 
the data suggest that further attention should be given 
to contaminants in VO in the future. On the other hand, 
the data indicate potential genotype-specific differences 
in the response of the intestinal transcriptome and 
proteome to dietary VO. These include potential 
changes in structural properties of the intestinal layer 
and defence against cellular stress suggesting the Lean 
group was more susceptible to diet-induced oxidative 
stress. Considering the possibility of selecting families 
better adapted to alternative diet formulations, and the 
central role of intestine as a major barrier to nutrients, 
contaminants and pathogens, greater attention should 
be given to this organ when evaluating the effects of diet 
and genotype. 

Methods 

Feeding trial and sampling 

A dietary trial was conducted using two genetically char- 
acterized groups of Atlantic salmon (Salmo salar) post- 
smolts comprising full-sib families selected from a 
breeding program (Landcatch Natural Selection Ltd., Ar- 
gyll, Scotland). The choice of the family groups was 
based on estimated breeding values (EBVs) of the par- 
ents for high or low flesh adiposity, assessed by Torry 
Fatmeter (Distell Industries, West Lothian, UK), a trait 
that was found to have a heritability ranging from 0.17 
to 0.39 in this dataset. The two groups were created 
from four unrelated full-sib families; two families from 
the extreme lower end of the EBV distribution for flesh 
lipid content (Lean) and two families from the extreme 
upper end of the distribution (Fat). The average EBV for 
the lipid content of the Fat families was 2.00 percentage 
units higher than that of the Lean families, representing 
a standardised selection differential of 2.33 standard 
deviations. Assessment of the flesh and visceral lipid 
contents at the end of the trial confirmed differences in 
adiposity between the groups [5]. 

Two thousand fish of each group were stocked into 
eight 12 x 5m 3 net pens at the Ardnish Fish Trials Unit 
(Marine Harvest Scotland, Lochailort, Scotland; 500 fish 
pen" 1 ). Duplicate pens of each group were fed one of 
two experimental diets (Skretting ARC, Stavanger, Nor- 
way) containing 25-32% fish meal, 40-45% plant meals 
and 27.5-30% oil (percentages varying according to pellet 
size) supplied either as standard northern FO or as a VO 
blend comprising rapeseed, palm and Camelina oils in a 
ratio of 5:3:2. Diets were formulated to fully satisfy the 
nutritional requirements of salmonid fish [79] and con- 
tained similar levels of PUFA but different n-3 and n-6 
PUFA contents, 25.3% and 4.6% in the FO diet and 
13.4% and 17.1% in the VO diet, respectively [5]. After 
55 weeks, 25 fish per pen were sampled 24 h after the 



last meal. Fish were killed by a blow to the head follow- 
ing anaesthesia, and intestinal tissue (pyloric caeca) col- 
lected, immediately frozen in liquid nitrogen and stored 
at -70°C prior to analyses. Further details can be found 
in Bell et al. [5]. 

Lipid extraction and fatty acid analyses 

Total lipid from 1 g of intestine of four fish per treat- 
ment was extracted and determined gravimetrically [80], 
and fatty acid methyl esters (FAME) prepared by acid- 
catalysed transesterification of total lipid [81]. FAME 
were separated and quantified by gas chromatography as 
described in detail previously [5]. Significant differences 
in intestinal fatty acid composition were determined by 
two-way ANOVA (p<0.05) using the SPSS 16.0 statistical 
package (SPSS Inc., Chicago IL, USA). 

RNA extraction and purification 

Intestinal tissue (0.2 g) from six individuals per experi- 
mental group was homogenised in 2mL TRI Reagent 
and total RNA isolated following manufacturer s instruc- 
tions (Ambion, Applied Biosystems, Warrington, U.K.). 
RNA quantity and quality (integrity and purity) were 
assessed by gel electrophoresis and spectrophotometry 
(NanoDrop ND-1000, Thermo Scientific, Wilmington, 
U.S.A.), and 100 ug of total RNA from each sample fur- 
ther cleaned by mini spin-column purification (RNeasy 
Mini Kit, Qiagen, Crawley, UK). 

Microarray hybridizations, image processing and 
statistical analysis 

The TRAITS/SGP (v.2.1) salmon 17k cDNA microarray 
described by Taggart et al. [33] was used (ArrayExpress 
accession: A-MEXP-1930). A dual-labelled experimental 
design was employed, with each sample being competi- 
tively hybridised against a common pooled-reference. 
The experiment comprised 2 genotypes (Lean/Fat) x 2 
diets (VO/FO) x 6 biological replicates. Indirect labelling 
was employed for preparing the microarray targets. 
Antisense amplified RNA (aRNA) was produced from 
500 ng of purified total RNA per sample using the 
Amino Allyl MessageAmpTM II aRNA Amplification 
Kit as per manufacturer s instructions (Ambion, Applied 
Biosystems), followed by Cy3 (for samples) or Cy5 (for 
pooled reference) fluor (GE HealthCare, Buckingham- 
shire, U.K.) incorporation through dye-coupling reaction 
[10]. Microarray hybridizations were performed in a 
Lucidea semi-automated system (GE Healthcare) with- 
out pre-hybridization. For each array, every labelled bio- 
logical replicate and corresponding pooled reference (40 
pmol each dye, c. 150 ng aRNA) were combined and added 
to the hybridization solution. Two post-hybridization 
automatic washes followed by six manual washes to a final 
stringency of 0.1 x SSC (EasyDipTM Slide staining system; 
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Canemco Inc., Quebec, Canada) were performed before 
scanning. 

Scanning was performed at 10 \xm resolution using an 
Axon GenePix 4200AL Scanner (MDS Analytical Tech- 
nologies, Wokingham, Berkshire, U.K.). Laser power was 
constant (80%) and "auto PMT" was enabled to adjust 
each channel at less than 0.1% feature saturation and 
Cy3/Cy5 mean intensity close to one. BlueFuse software 
(BlueGnome, Cambridge, U.K.) was used to identify fea- 
tures and extract fluorescence intensity values from TIF 
images. The resulting fluorescence intensity data (Blue- 
Fuse proprietary algorithm) and quality annotations for 
the 17,102 gene features were exported into the Gene- 
Spring GX version 10.0.2 analysis platform (Agilent 
Technologies, Wokingham, Berkshire, U.K.) after under- 
going block Lowess normalization. All control features 
were excluded from subsequent analyses. Data trans- 
formation and quality filtering were as in Morais et al. 
[9,10]. This gave a final list of 15,498 genes that were eli- 
gible for statistical analysis. Experimental annotations 
complied fully with minimum information about a 
microarray experiment (MIAME) guidelines [82] and ex- 
perimental hybridisations are archived on the EBI 
ArrayExpress database (http://www.ebi.ac.uk/arrayexpress/) 
under accession number E-TABM-1173. Hybridization 
data were analysed in GeneSpring by two-way ANOVA, 
which examined the explanatory power of the variable 
diet and genotype and the interaction between the two, 
followed by Gene Ontology (GO) enrichment analysis of 
the significant lists of features, at a significance level of 
0.05. No multiple test correction was employed, as pre- 
vious analyses, confirmed by RT-qPCR, indicated that 
such corrections are over-conservative for this type of 
data [9,32]. 

RT-qPCR gene expression analysis 

Expression of selected genes, for microarray validation 
and to further examine biological processes of interest, 
was studied by reverse transcription quantitative real 
time PCR (RT-qPCR), with target qPCR primer 
sequences given in Additional file 2. In addition, amplifi- 
cation of two reference genes, cofilin-2 and elongation 
factor- la {elf- la), was performed. One \ig of column- 
purified total RNA per sample was reverse transcribed 
into cDNA using the Verso TM cDNA kit (ABgene, Sur- 
rey, U.K.) using a mixture of random hexamers (400ng/ 
\iL) and anchored oligo-dT (500ng/|iL) at 3:1 (v/v). 
Negative controls (no enzyme) were performed to check 
for genomic DNA contamination. A similar amount of 
cDNA was pooled from all samples and the remaining 
cDNA diluted 20-fold with water. RT-qPCR analysis used 
relative quantification with the amplification efficiency of 
each primer pair assessed by serial dilutions of the cDNA 
pool. Amplifications were carried out in duplicate using 



a Quantica machine (Techne, Cambridge, U.K.) in a final 
volume of 20 [A containing 2-8 [A diluted cDNA, 0.5 [iM 
of each primer and 10 [A AbsoluteTM QPCR SYBR® 
Green mix (ABgene), with a systematic negative control 
(NTC-non template control). The qPCR profiles con- 
tained an initial activation step at 95°C for 15 min, fol- 
lowed by 30-40 cycles: 15 s at 95°C, 15 s at the specific 
primer pair annealing temperature (Ta; Additional file 2) 
and 15 s at 72°C. After amplification, a melt curve was 
performed confirming a single product in each reaction, 
RT-qPCR product sizes checked by agarose gel electro- 
phoresis, and identity of amplicons confirmed by sequen- 
cing. Gene expression was analysed using the relative 
expression software tool (REST 2008, http://www.gene- 
quantification.info/), employing a pair wise fixed reallo- 
cation randomisation test (10,000 randomisations) with 
efficiency correction [83]. 

Protein extraction and labelling 

Six intestine samples (100 mg) per treatment were rap- 
idly disrupted by homogenization and sonication on ice 
in 1 ml of DIGE lysis/labeling buffer (7M Urea, 2M Thio- 
urea, 4% CHAPS and 30mM Tris, pH 8.5) in the pres- 
ence of 10 \A of a protease inhibitor cocktail (ref. P8340, 
Sigma- Aldrich) and 4[il of 250 mM EDTA. After centri- 
fugation at 12,000 x g for 20 min at 4°C, the supernatant 
was recovered and protein concentration determined 
(Quick Start Bradford Protein Assay Kit, ref. # 500-0202, 
Bio-Rad). Protein (500 \ig) was purified by precipitation 
(ReadyPrepTM 2-D Cleanup kit, Bio-Rad) and the pellet 
re-suspended in DIGE lysis/labeling buffer at ~5|ig/|il. 

Samples were labelled using CyDye DIGE fluors 
(5nmol labelling kit; GE Healthcare), following manufac- 
turer s instructions. Three of the experimental replicates 
of each treatment (50 \ig protein with pH corrected to 
8.0-8.5) were labelled individually with 400 pmol Cy3 
and the remaining three with 400 pmol Cy5. In addition, 
equal amounts of all experimental samples were pooled 
and 600 \ig of protein (50 \xg per 2-D gel) were batch- 
labelled with Cy2. The three labelled samples, corre- 
sponding to two experimental samples and one internal 
reference pool, were then combined to have in each 2-D 
gel samples corresponding to fish fed either FO or VO 
within the same family group. 

Two-dimensional (2-D) polyacrylamide gel electrophoresis 

Rehydration buffer (ReadyPrep 2-D starter kit, Bio-Rad) 
containing 0.2% DTT was added to the pooled protein 
samples to a final volume of 450 [A, which were loaded 
onto Immobiline DryStrip pH 3-11 NL, 24 cm IPG strips 
(GE Healthcare) by passive rehydration at room 
temperature overnight in the dark. Proteins were sepa- 
rated in the first dimension by isoelectric focusing at 20°C 
(Ettan IPGphor, GE Healthcare), applying increasing 
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voltage until 200 V for 4 h, increasing to 500 V over a 
period of 3 h, then keeping the applied tension at a con- 
stant 1000 V for 1 h, followed by a further increase to 
8000 V over 90 min, maintaining this voltage for almost 
9 h (total of 79,200 Vh). After isoelectric focusing the 
strips were equilibrated in two 40 min steps using 
50mM Tris-HCl pH 8.8, 6M urea, 30% glycerol, 2% SDS 
buffer, to which 2 % DTT (w/v) and 2.8% iodoacetamide 
(GE Healthcare) were added to produce reducing and al- 
kylating buffers, respectively. The strips were loaded 
onto a 12.5% acrylamide gel cast between low fluores- 
cence glass cassettes (EttanDALTsix gel caster system, 
GE Healthcare). The strips were overlaid with ReadyPrep 
Overlay Agarose (Bio-Rad) and the six gel cassettes run 
in the EttanDALT system in two steps: at 60 mA, 80 V, 
6 W (upper limits) for 1 h, and then 240 mA, 500 V, 
78 W until the bromophenol blue dye front had run to 
~1 cm above the bottom of the gels. Laemmli buffers 
(250mM Tris, 1.92M glycine, 1% SDS) were used in the 
lower (IX) and upper (2X) chambers, respectively. 

Gel imaging and analysis 

Labelled gels were scanned using a Typhoon TRIO (GE 
Healthcare) and Cy2, Cy3 and Cy5 images acquired 
using 520BP40, 580BP30 and 670BP30 laser emission fil- 
ters, respectively, at 500 PMT and 100 um resolution. 
Images were cropped to remove extraneous areas prior 
to analysis (ImageQuantTL, GE Healthcare), and image 
analysis performed using DeCyder V7.0 (GE Healthcare). 
The estimated number of spots for each co-detection 
procedure was set at 10,000 and an exclusion filter was 
applied to remove spots with a volume lower than 
30,000. Differential expression of protein spots was 
examined by two-way ANOVA at a significance level of 
0.05. After verifying that significant spots were well 
matched across the gels, two pick lists were generated 
with a total of 22 and 45 spots for the diet and genotype 
factors, respectively (Figure 1; Additional files 3 and 4). 

Spot picking and protein identification by peptide 
fragment fingerprinting 

Four preparative gels were run under the conditions 
described above but with higher amounts of protein (600 
and 1000 ug). They were stained with colloidal Coomas- 
sie and, whenever possible, spots were excised and 
sequenced in the Mass Spectrometry Laboratory ITQB- 
UNL (Oeiras, Portugal), where in-gel digestion and ex- 
traction of the proteins from the gel was performed, fol- 
lowed by micropurification, and peptides identified by 
mass spectrometry (using an Applied Biosystems (ABI) 
4800 MALDI TOF/TOF Analyzer). The search engine 
MASCOT was then used to identify and confirm protein 
ID s from the peptide mass fingerprinting (PMF; MS) and 
peptide fragment fingerprinting (PFF; MS/MS) data. 



Additional files 



Additional file 1: Intestine transcripts corresponding to the top 100 
most significant features exhibiting a significant diet x family 
interaction. 

Additional file 2: Primers used for RT-qPCR analyses [84 89]. 
Additional file 3: 2D-gels pick list for the factor diet. 
Additional file 4: 2D-gels pick list for the factor genotype. 
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